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Abstract 

A "toy" model, simple and elementary enough for an undergraduate class, of the temperature 
dependence of the greenhouse (mid-IR) absorption by atmospheric water vapor implies a bistable 
climate system. The stable states are glaciation and warm interglacials, while intermediate states 
are unstable. This is in qualitative accord with the paleoclimatic data. The present climate may 
be unstable, with or without anthropogenic interventions such as CO2 emission, unless there is 
additional stabilizing feedback such as "geoengineering" . 
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I. INTRODUCTION 



Models of climate change, a subject that concerns time scales from decades to billions 
of years, range from the simple and qualitative to atmospheric and oceanic general circu- 
lation models running on the fastest computers. Despite several decades of effort by many 
scientists, the causes of such important phenomena as the alternation of glaciation and in- 
terglacials, and the ice ages themselves, remain controversial. In this pedagogical paper I 
describe an elementary and qualitative model of the consequences of one important effect, 
the temperature dependence of the opacity of atmospheric water vapor, the most important 
greenhouse gas 1 . This is a consequence of the temperature dependence of the water vapor 
pressure and the physics of pressure broadening of saturated Lorentzian spectral line profiles, 
and lies on the border between physics and geophysics. 

II. PALEOCLIMATIC DATA 

The paleoclimatic data are summarized in Figure [TJ Alternation of ice ages and inter- 
glacials is evident. Intermediate states, such as the present climate, are not stable or steady, 
but show continual variation. This is also evident in the historical record of the Late Me- 
dieval Climatic Maximum (c. 800-1300), followed by the Little Ice Age (c. 1300-1800), 
followed in turn by a warming trend. 

III. WATER VAPOR OPACITY 

Although water vapor is the greatest source of atmospheric opacity in the mid-infrared (8- 
14 n) band in which the Earth's thermal emission peaks, it is not even mentioned in a classic 
monograph 6 on climatological dynamics. The HITRAN database^ is the comprehensive 
source of spectral information. Measurements of atmospheric opacity, such as those shown 
in Fig. 3 of Hiriart and Salas^, illustrate its sensitivity to even small amounts of water vapor. 

Tropospheric water vapor is not usually considered in the context of anthropogenic green- 
house gas (GHG) emissions because it equilibrates with surface water very rapidly, with a 
characteristic time of a few days or less, and any anthropogenic contribution has no long 
term effect. Tropospheric water vapor opacity is controlled by evaporation and precipitation 
to a level determined by the surface temperature. 
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Ice Age Temperature Changes 




IV. LIMIT STATES 

This sensitivity of opacity to water vapor pressure, combined with the sensitivity of 
equilibrium vapor pressures to temperature, qualitatively suggests a bistable system: There 
may be a cold state in which the water precipitates and the atmosphere is transparent, with 
little water vapor absorption, and a warm state in which the atmosphere is warm, humid and 
opaque and in which absorbed solar energy is carried through the troposphere by convection 
(resembling present tropical or temperate summer conditions). 

These limit states are stable: in the cold state perturbations to the temperature have 
little effect on the infrared opacity because it remains very low, while in the warm state the 
opacity has little effect on the heat transfer. In either case there is little feedback to the heat 
balance from changes in the tropospheric water vapor content. At intermediate temperatures 
increases in temperature produce a positive feedback by increasing the infrared opacity (and 
vice versa for decreases), so that perturbations grow exponentially until one of the stable 
limit states is approached. 
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V. MODEL 



Quantitative calculation is precluded by the complexity of the climate system. These 
complexities includ processes that are difficult to calculate but understood in principle (such 
as opacity resulting from a dense comb of infrared vibration-rotation lines) and those not 
understood (such as cloud formation and greenhouse gas exchange with the ocean, soil and 
biosphere). Even the most elaborate and sophisticated general circulation models (GCM) 
contain many uncertain parameters that can only be calibrated at a single measured point: 
today's climate. It is not surprising that different GCM disagree in the magnitude of the 
effect of GHG emission, and it is not evident how to resolve these disagreements. 

Disagreements among the best state-of-the-art GCM suggest considering the opposite 
approach, a model naive almost to the point of triviality. If its workings are transparent, it 
may provide useful insight, even though it cannot make quantitative predictions. 

Consider a surface heat reservoir warmed by sunlight with a mean (diurnally and annually 
averaged) intensity P Q and cooled by black body radiative emission with a (Planck weighted) 
fraction fh(T) of the infrared spectrum blocked by lines of Lorentzian lineshape^ (Doppler 
broadening of molecular lines is negligible at atmospheric temperatures). The optical depth 
at the centers of strong lines tq 3> 1, and assume an underlying heat capacity (per unit area) 
C p . The temperature T of the thermal reservoir changes at a rate 



where <tsb is the Stefan-Boltzmann constant and we approximate the line profile as com- 
pletely opaque within a blocked fraction fb(T) of the spectrum, and completely transparent 
outside it. 

We suppose an initial equilibrium temperature Tq defined by the heat balance condition: 




(1) 




(2) 



Small deviations of temperature are described by 



T = T + ST. 



(3) 
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VI. INSTABILITY 



Expanding Eq. [T] to first order in ST: 

1 dST 



u T dt 



ST 



where 



vsbT[ 



ST grows or decays with an exponentiation rate v: 

v ~ dh(T) 
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vq dT 

To estimate the derivative in Eq. [HI we use an approximate equation^ for the saturation 
water vapor pressure over a surface of temperature T: 



V: 



,(T) oc T 3/2 exp (-AH/k B T), 



(7) 



where AH = 2260 J/g is the enthalpy of evaporation of water and k B is Boltzmann's 
constant. For a Lorentzian line profile with To 3> 1, the width and ft, at a specified optical 
depth 0(1) 
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where we have made the tacit assumptions 1 — fb(T ) > r , which amount to assuming 
that there is enough unblocked spectrum for the line to broaden according to Eq. |HJ The 
appropriate value of tq in this condition is an average over all lines that are optically thick at 
their cores. This is likely to be dominated by the very numerous lines for which r exceeds 
unity by a factor of a few, but not by orders of magnitude, so these results may be valid for 
fb ^ 0.7. The model breaks down as /& — >■ 1 because then there is little unblocked spectrum 
into which the lines can broaden. 

Evaluating at a sea surface T = 290 °K, 



- = 13.2/ 6 (T )-4. 



(11) 



A quantitative calculation of /& is not possible within this qualitative model, but spectral 
data 8 suggest values of several tenths. Then deviations from equilibrium grow exponentially 
for 



The mean ocean depth, averaged over the entire Earth (including zero depth for land 
areas) is about 2.65 km. However, the flow of heat to the deep ocean, while apparently 
reflecting the warming of the last century 11 , is expected to be confined to downwelling 
regions and to be comparatively slow. For the purpose of an order of magnitude estimate 
we suppose a nominal mean depth of 1 km contributes to C p and find 



VII. CONCLUSION 

If 0.3 < fb < 0.7, as suggested by the empirical water vapor spectrum 8 -, then instabil- 
ity would be expected for any state between an ice age (when p v is low enough that the 
assumption tq ^> 1 is not satisfied) and the warmest interglacials (when the assumption 
1 — /b(To) > t 1 ' 2 is not satisfied as nearly the entire thermal infrared spectrum is blan- 
keted with water lines and /& — > 1, or energy is carried by convection, ignored in this naive 
model). Hence the climate would be stable only in these two limiting states. 

The paleoclimatic data, as shown in Figured! are qualitatively consistent with this model. 
The historic record of the Medieval Climatic Maximum, the Little Ice Age, and fluctuations 
on similar time scales of (9(100) years at other epochs, is also qualitatively consistent with 
the estimated u Q . 

VIII. DISCUSSION 

Climatologists have struggled with the question of what upsets the stable ice age and 
interglacial states since the discovery of the ice ages, and this naive model cannot contribute 
to its answer. Even the most complete and sophisticated general circulation models struggle 
with this problem, and its resolution remains controversial. However, even the conclusion 
that intermediate states may be intrinsically unstable has significant implications. One 
implication is that there is no "normal" climate, nor can the Earth be expected to return 



0.3 <f b < 0.7. 



(12) 



u Q = 1.04 x 10" 2 /y. 



(13) 
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to such a state even without anthropogenic interventions such as the emission of GHG. A 
second implication is that if some climate state is determined to be optimal for humanity, 
maintaining it would require continual active intervention (geoengineering) . 
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